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The particles of a dark matter due to gravitational interaction deviate from straight trajectories in the 
vicinity of a massive body. This causes their density to become inhomogeneous. The developed density 
contrast causes a gravitation attraction force acting upon the body proportional to its mass squared. Since 
this effect depends on the body's velocity, it causes a specific reference frame to stand out. The force is 
similar to an anisotropic drag, which can be negative in some directions. The mean drag force, averaged 
over all direction, is nonpositive. We can expect some observational manifestation of the considered effect for 
supermassive black holes in galaxies. 
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I. INTRODUCTION 

The last couple of decades essentially changed our con- 
cept about the Universe and its contents. The ordinary 
baryonic matter amounts to 4.6% of its mean density, 
72.6% corresponds to the dark energy and 22.8% to the 
dark matter according to the results of data processing 
of Wilkinson Microwave Anisotropy Probe observations 
of the cosmic microwave background (CMB) radiation 
anisotropy Astronomers predicted the existence of 
dark or unseen matter starting from the 1933 Zwicky pa- 
per Now we have a lot of evidence from astronomy, 
astrophysics, and cosmology in support of this fact. Let 
us indicate some of them. A large difference between 
the virial masses of galaxy clusters and of the total mass 
of all galaxies belonging to these clusters indicates the 
presence of dark matter in clusters. Rotation curves of 
galaxies show that the major part of their masses belongs 
to a dark halo rather than to a stellar disk. We can also 
mention merging galaxy clusters, e.g. the Bullet Clus- 
ter 1E0657-558, strong and weak gravitational lensing, 
baryonic acoustic oscillations, simulations of galaxy and 
large-scale structure formation, large-scale collective mo- 
tion of galaxies, and type la supernovae explosions. The 
most important, however, is the CMB anisotropy data, 
which indicates the nonbaryonic nature of the majority 
of dark matter. 

There are lots of dark matter candidates such as ster- 
ile neutrinos, axions, neutralinos, gravitinos, and weakly 
interacting massive particles, which are explained in the 
review Q. The dark matter must interact gravitation- 
ally and cannot take part in electromagnetic or strong 
interactions. It is possible that the dark matter can 
also experience weak interaction. The experiments for 
the dark matter direct search such as Xenon 10/100/1T, 
DAMA/LIBRA, CDMS, CoGeNT, CRESST and many 
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others are trying to detect the dark matter assuming that 
it interacts weakly with ordinary matter. These experi- 
ments are explained in the article [4| . They constrain the 
cross section of weak interaction of dark matter to very 
small values. Nevertheless, there is a possibility that the 
dark matter cannot interact with the ordinary matter via 
the weak force at all. Sometimes the dark matter, which 
cannot interact weakly is referred to as a mirror matter. 

We consider only the gravitational interaction between 
the particles of the dark matter with bodies made from 
the ordinary matter or black holes. We exclude the weak 
interaction since we are only interested in the effect due 
to gravity. Naturally, there is lot of known effects of the 
gravitational attraction between such objects and mas- 
sive aggregates of dark matter. We do not consider them. 
Instead, we concentrate on the gravitational interaction 
with the dark matter particles flying through and near 
an object. The considered effect arises due to the trajec- 
tory bending of dark matter particles in the gravitational 
field of the object and the gravitational attraction of the 
object to these particles. 

We estimate the value of the force and the acceleration 
caused by this effect and find the class of the objects 
for which this effect is essential. We consider the cases 
of isotropic and anisotropic velocity distributions of the 
dark matter particles and show that in both cases there 
is a specific reference frame, which stands out due to this 
effect. 



II. CAUSE OF THE EFFECT 

Let us start from a preliminary problem. A spherical 
body (e.g. the Sun or any other star) has the mass M 
and the radius R. It rests at the origin of the cylindrical 
coordinate system (p,4>,z). A flow of the dark matter 
passes through and around the body. Far from the body 
at z <C —R the flow is homogeneous and moves along the 
z axis with the initial velocity v. Our goal is to calculate 
the force acting on the body due to the gravitational 



attraction to dark matter particles. 

We consider the case of cold dark matter with the ve- 
locity much less than the speed of light v <C c and use 
the Newtonian mechanics and Newton's law of univer- 
sal gravitation. The masses of dark matter particles are 
negligibly small as compared to M. 

Particles' trajectories deviate from straight lines in the 
gravitational field of the central body. The majority 
misses the body and moves along hyperbolae. The tra- 
jectory of a particle deviates by the angle ip 



ip = 2 arctan 



MG 
v 2 p ' 



(1) 



Here po > R is an impact parameter and G is the grav- 
itational constant. We assume this angle to be small 
ip <C 1, so we can assume ip ~ 2MG/(v 2 po). In this case 
the distance from the center of the body to a particle's 
periapsis or a point of closest approach is almost equal to 
the impact parameter po- As it will be shown further, the 
particles with p > R provide a lion share of the effect. 
The particles flying through the body with po < R make 
a small contribution to it, so we can restrict ourselves 
to a rough estimation. The deviation angle ip for such 
particles depends on the density distribution of baryonic 
matter inside the body. We know that at po = we have 
ip = from the axial symmetry and at po — R this angle 
must match So we consider an approximation 
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(2) 



with some positive exponent 7 > 0. 

In order to demonstrate the physical nature of the ef- 
fect we consider at first a simple toy model (see Fig. 
[T]). We substitute the real hyperbolic trajectories of 
dark matter particles (the dashed curve in Fig. [1]) by 
their asymptotes. So in the framework of this toy model 
particles move along straight lines p = po parallel to 
the z axis at z < 0, at z — they turn by the an- 
gles V'(Po) an d continue moving along the straight lines 
p = \po — vt sin?/; I, z = vt cos-0, where t — corresponds 
to the moment of crossing the plane z — 0. These lines 
intersect at the axis p = and then diverge. In addition 
we consider the motion with constant velocity v. 

The particle density in the half-space z < is constant. 
The mass dm passing through the area dS of the plane 
z = const < during the time interval dt is 



dm = AdSdt. 



(3) 



Here A is the flux of the dark matter. The penetrat- 
ing particles occupy the volume dV = dSvdt, thus the 
dark matter density pdm can be expressed as pdm = 
dm/dV = A/v. In the half-space z > the distribu- 
tion of particle density becomes inhomogeneous due to 
the trajectories bending. This leads to the considered 
effect, namely the difference between the gravitational 
attractions between the body and the particles in both 




FIG. 1. Dark matter particles in the gravitational field of the 
spherical body 



half-spaces. To calculate the resultant force let us con- 
sider the particles passing through a ring-shaped area 
with po < p < po + dp during the time interval from 
t = to t = dt. At time t > these particles form a ring 
with the center at the point p = 0, z = vtcosip. The dis- 
tance between the particles and the center of the body is 

1/2 

r = (p 2 — 2povt sintp + v 2 t 2 ) . The mass of the parti- 
cles according to ([3]) is dm — ■nAp^dpodt. The force due 
to the gravitational attraction between the body and the 
particles is aligned with the z axis and equals to 



GMzdm 



= 2-kpdmGMv 2 I cosippodpo 



tdt 



3/2 



(4) 



(p 2 — 2povt sintp + v 2 t 2 ) 



= 2irp DM GM I 1 + 8 ™ dp 
cost/; 



The force of the gravitational attraction between the 
body and dark matter in the half-space z < can be 
obtained in a simple way. It is equal to (HJ in the case 
when the integral is taken between the time limits — 00 
and at ip = 0. Thus its z component is equal to 
F_ = -F + (tp = 0). The resultant force F = F + + F_. 
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Using Q we get for the case of small deviation angles ip 



F = 2-kpdmGM / ipdp 

R 
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This integral diverges at large po values, thus we intro- 
duced a finite upper limit R m ax- 

The same result can be also obtained straightforwardly 
without any simplification of the calculations. The dark 
matter particles missing the body move along hyperbolae 
with varying velocities according to the laws of orbital 
mechanics. We can introduce the useful parameterization 
t = GMv~ 3 {e cosh£ - £), r = GMv~ 2 {e cosh£ - 1), z 
r sin(/0/2 + arccos /^g^g'-i ) w ^h deviation angle ip from 
([1]) and the eccentricity of the hyperbola e = 1/ sin(^/2). 
After integration with respect to £ between the limits 
— oo and oo we obtain for the case C 1 exactly the 
expression ([5]). 

Let us note some similarity of this force and the ion- 
ization losses of highly energetic charged particles or the 
dynamical friction due to the gravitational interaction of 
a massive, rapidly moving star with a cluster of stars 
(see In all cases F oc v~ 2 \n(R max / R min ) and the 
formula for the dynamical friction coincides with ([S]). 
Nevertheless, the considered problem differs from the dy- 
namical friction. In the case of the dynamical friction we 
consider a fast body gravitationally interacting with slow 
ones; in our case we consider a slow or immobile body 
gravitationally interacting with fast particles. If these 
particles move with the same initial velocity, we can use 
the frame of reference moving with this initial velocity 
and obtain the situation similar to the dynamical fric- 
tion case. Naturally, the force must be the same and this 
condition is held. But in the general case in which parti- 
cles move in different directions with different velocities 
we cannot reduce the considered problem to the dynam- 
ical friction. Later on we will consider this more realistic 
general case and find that the energy of the body can 
not only decrease, but also increase due to the interac- 
tion with particles. 



III. ESTIMATION OF THE VALUE OF THE FORCE 

The force F is proportional to the square of the mass 
of the body, so it provides an acceleration proportional 
to the mass 



A-k Pdm G 2 M R 

max- 
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R 



(6) 



This acceleration is very small because it is propor- 
tional to the product of two small parameters, namely 
the gravitational constant squared and the dark matter 
density. It is negligible for ordinary bodies, therefore this 
effect cannot be found in direct experiment. If we draw 
attention to compact astronomical objects, we shall find 
a lot of bodies with big masses. Moreover, the astro- 
nomical time scale provides a lot of time for small but 
permanent acceleration to manifest. 

Let us estimate by an order of the magnitude the 
acceleration of the Sun according to ©• Its mass is 
M S un ~ 2 x 10 30 kg and the radius R Sun « 7 x 10 s m. 
We can estimate R ma x as the half-distance to the nearest 
stars Rmax ~ 1 ly ~ 10 16 m thus \n(R max /R) pa 16. Ac- 
cording to the modern paradigm, the dark matter forms a 
dark halo around the Milky Way galaxy and provides the 
majority of its mass pa l0 12 Mg un . Taking the diameter of 
the dark halo as 500,000 ly pa 5 x 10 21 m we can estimate 
Pdm ~ 2 x 10~ 22 kg m~ 3 . This is an underestimated 
value because of the inhomogeneity of the dark matter 
distribution. Using the Navarro-Frenk- White fitting for- 
mula [6( for it we can improve this rough estimation. Also 
we can use the estimations based on astronomical obser- 
vations. A number of pdm estimations are cited in the 
review They vary in the interval 0.2-^0.8 GeV cm -3 . 
This range was confirmed later by Weber and de Boer 
0- Choosing the value p DM ~ 0.5 GeV cm -3 pa 10~ 21 
kg m~ 3 and assuming v pa 200 km s _1 we obtain the 
estimation a pa 2 x 10~ 20 m s~ 2 . 

This acceleration is less than the constant ao in the 
MOND theory by 10 orders of magnitude. During the 
period of the rotation of the Sun around the galactic 
center pa 250 Myr ss 8 x 10 15 s this acceleration changes 
its velocity only by 0.16 mm s _1 . During the age of the 
Universe 14 Gyr it changes the velocity of the Sun by 0.9 
cm s _1 . 

Thus we see that there is a theoretical possibility to 
transfer the energy, the momentum and the angular mo- 
mentum of the mirror dark matter to the ordinary mat- 
ter. Nevertheless this effect is very small and virtually 
cannot affect the motion of usual celestial bodies. 

Could this effect become significant for other bodies 
or under different circumstances? We used the minimal 
possible assumption of the value of the typical dark mat- 
ter particles velocity. It corresponds to the Sun's orbital 
velocity around the center of the Galaxy. Thus we cannot 
decrease the denominator in ([6]). We have two possibili- 
ties to increase the numerator. The dark matter density 
was essentially greater in the early Universe, but there 
were no compact bodies. Thus there remains only one 
thing to consider - objects with the mass much greater 
than Msun- Stars have masses up to ~ 100 Ms-am but 
the most massive ones have lifetimes about a few million 
years. So the effect is small also for the most massive 
stars. If we consider the acceleration of galaxies we have 
to use for pnu the cosmological value pa 0.2p c pa 2 x 10~ 27 
kg m~ 3 , where p c is the critical cosmological density of 
the Friedmann universe. Thus, the acceleration is small 



for galaxies as well. 

The most promising objects are supermassive black 
holes, e.g. the Sagittarius A* black hole with the mass 
4 x 10 6 Msun in the center of our Galaxy Q. Even with- 
out taking into account the greater dark matter density 
in the galactic centre we estimate a « 8 x 10~ 14 m s~ 2 . 
During the period of the rotation of the Sun around the 
galactic center its velocity would change by 640 m s _1 . 
Black holes in galaxies, e.g. active galactic nuclei (AGN) 
can be much more massive. The largest supermassive 
black hole in the Milky Way's neighborhood appears to 
be that of M87 galaxy, weighing about 6 x 10 9 Ms un @ 
and the most massive known black hole in the NGC 4889 
galaxy has the mass 2 x 10 10 Ms U n [13] ■ During the pe- 
riod of the rotation of the Sun around the galactic center 
the velocity of the latter changes by 3200 km s _1 . So the 
considered effect can in principle have some astronomical 
manifestations, but only for supermassive black holes. 

There is another effect, specific for black holes, which 
arises due to the capture of dark matter particles by black 
holes. We will demonstrate that it is much weaker than 
the above mentioned one below. 



IV. ISOTROPIC CASE 




FIG. 2. Relations between velocities and angles in the rest 
and in the moving frames 



Let us return to the formal problem. There is no rea- 
son to assume that the dark matter comes only from one 
direction. We should have used rather a different prob- 
lem formulation even expecting one preferential direction 
of dark matter motion. In the one-dimensional version of 
the problem the dark matter approaches the body along 
the direction of the z axis. We suppose that there is a 
frame of reference in which the velocities and the densi- 
ties of the particles approaching from two opposite direc- 
tions are equal. We denote the particles velocity by u. 
The body is moving along the z axis with the velocity vo, 
\vq\ < u. The gravitational force due to the interaction 
with dark matter moving along the z axis is directed to- 
wards it and characterized by ([5]) with v = u — vq. The 
gravitational force due to the interaction with dark mat- 
ter moving along the z axis in the negative direction is 
directed opposite to the first force and characterized by 
([5]) with v = u + vq. The resultant force in the direction 
of the axis is 

F = ±4irp DM G 2 M 2 In 

o /~<2 a ,r2 l Umax 

= 8irp DM & M In — — , , , 
R (u z — vo ) 

The factor 1/2 arises because the dark matter density 
duplicates for the same A due to the two possible direc- 
tions of approaching. For a small velocity of the body this 
force is proportional to vq similar to the drag, but with 
negative viscosity. The value of vq increases with time. 
At small velocities |uq| <C m it increases exponentially. 

Does this effect survive if we consider more realistic 
models in which dark matter particles approach from all 
directions? Let us begin with the simplest one. Let us 
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assume that there is a reference frame in which the dark 
matter flow is isotropic. The particles in this frame have 
the velocity u and their directions are distributed uni- 
formly. Introducing a spherical coordinate system with 
angular coordinates 9 and </>, we can consider the mass of 
the dark matter dm passing through the area dS of some 
surface during the time interval dt from the spatial angle 
dQ in the form similar to the definition of luminance in 
optics, 



dm = B( 



cos (3 dS dt dfl. 



(8) 



Here /3 is the angle between the particles' direction of 
arrival and the direction of the surface normal. In 
the isotropic case i?=const. It can be easily expressed 
through the dark matter density. The density of the par- 
ticles with velocity directions in the spatial angle d£l is 
dp dm = Bu~ 1 dfl. After integration with respect to dfl 
we obtain pdm = ^irBu^ 1 . 

Particles moving with the velocity u from the spa- 
tial angle dVl provide the gravitational force affecting the 
body in this frame, which we will call the "rest" frame, 



dF= kd£DM^ k = 4nG 2 MHn R 



u 

u 



u 

dpDM 



BdQ 



R 



(9) 



Naturally, after the integration with respect to all direc- 
tions of the particles arrival we obtain F = 0. 

Let us consider the body moving with the velocity 
v with respect to the rest frame with the assumption 
vq < u. We introduce the ratio A = vq/u < 1. We will 
use the "moving" reference frame in which the body is 
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motionless. The velocities of the particles in this system 
v and their components are linked with the ones in the 
rest system u by simple relations [ll[ 



+ v , v = u + v + 2uv cost 
v cos a — vq + u cos 9, 
v sin a = u sin (9. 



(10) 



The orientation of Cartesian coordinates x, y is shown in 
Fig. [2j In the moving frame we introduced the spherical 
coordinates a and (f>. The angle (j) is the same as in the 
rest frame and the relation between the angles 9 and a 
follows from (fTOl). 



cos a 



A + cos ( 



(l + A 2 + 2Acos6») 1 /2^ 



cos 9 = cos a (1 — A~ sin 



>V2 



A sin 2 a. 



(11) 



We need the relation between the spatial angle in which 
the particles come from in the rest frame dfi and the 
one the same particles come from in the moving frame, 
which we denote d£l' (hereafter all angles with primes 
correspond to the moving frame) . Since d£l' = sin adad(j) 
and dil = sm9d9d<fi, so dtt — f(a)dil' with 



f(a) — 2A cos a + 



1 + A 2 cos 2a 
(1 - A 2 sin 2 ^ 1 ^' 



(12) 



The particles approaching to the body in the moving 
frame with the velocity v from the spatial angle dfl' pro- 
vide the gravitational force affecting the body 



dF 



kdp D M , B'(a,4>)dtt' 

^ — v,dp DM = • (13) 



The value dp dm must be the same in any frame of refer- 
ence. Comparing (ITU1) and (fT"3")) we obtain 



B{0,4>)v 
uf(a) 



(14) 



Let us calculate the force acting on the body in the 
moving frame in the case when the dark matter distribu- 
tion in the rest frame is isotropic B(9, (/>)=const. From 
the axial symmetry it is obvious that this force must be 
directed along the x axis. So we can restrict ourselves to 
calculating only F x , 



F, 



kup 



DM 



2 

kpDM 



cos a sin t 



2 i/ 



v{9) 2 

(A + cos 6) sin 9d9 
(l + A 2 + 2Acos(9) 3 /' 



(15) 



= 0. 



Thus in the isotropic case with vq < u the considered 
force is absent not only in the rest frame but also in any 
moving frame. 

We can use an analogy with the electrostatics to ex- 
plain this result. Let us consider the model similar to Fig. 
[21 but in the coordinate space, not in the velocity space. 



The thin spherical shell with radius r is charged with 
the distribution of the surface charge density <j{9,4>) pro- 
portional to B(6,4>). If we calculate the Coulomb force 
Fc acting on the test particle located in the point at 
x = — Ar, y = 0, then it will be proportional to the force 
F due to considered effect. The reason is that (p~3|) is 
similar to the Coulomb's law, but in the velocity space. 
In the isotropic case a = const and we deal with a uni- 
formly charged spherical shell. There is no electric field 
inside such a shell, so there is no force acting on the body 
in the moving frame. 

We can carry on this analogy. Instead of the case when 
the approaching particles have fixed velocities u in the 
rest frame, let us consider the case where the velocity 
has some distribution. The mass of the dark matter dm 
passing through the area dS of some surface during the 
time interval dt from the spatial angle dfl with the initial 
velocity v from the interval u < v < u + du is 



dm = B„. 



, u) cos ftdSdtdQdu. 



(16) 



In the analogy with the electrostatics u 3 B u (9,4>,u) play 
the role similar to the charge density p q . Instead of the 
thin spherical shell we deal with the charged ball. Its 
charge density p q (6, (f), r) is proportional to u 3 B u (9, <ft, u) 
and r is proportional to u. Let us consider the isotropic 
case B u (9,(j),u) = B u (u). Let the velocities of the ap- 
proaching particles in the rest frame be distributed in 
the range from u m i n to Umax- In the analogy with the 
electrostatics we deal with spherically symmetric charged 
thick shell with the inner radius r m i n cx u m i n and the 
outer radius r max oc u max . The case vq < u m i n corre- 
sponds to the configuration in which the test particle is 
located inside the cavity without the electric field. Thus 
there is no effect in any moving frame with vq < u m -m 
if the motion of the dark matter particles is isotropic in 
the rest frame. 

This conclusion completely changes if vq > u m i n . The 
force 



J B u (u)u 1 rfu 

F = V 

Vq 3 



(17) 



acts on the body in the "moving" frame. It corresponds 
to a negative friction. This effect takes place only at v > 
Umin- If u rn i n = this condition is satisfied automatically 
for any values. 



V. ANISOTROPIC CASE 

Let us consider the general case when the angular dis- 
tribution of arriving particles is anisotropic in any ref- 
erence frame. We start from the arbitrary initial ref- 
erence frame. Once again we utilize the electrostatic 
analogy and consider the continuous distribution of p q oc 
v?B u {9,<f),u) > over the velocity space. This charge 
distribution, which we assume to cover a finite area in the 
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velocity space, produces some electric field. The distri- 
bution of the field intensity over the velocity space must 
have at least one singular point at which its intensity is 
equal to zero. The existence of such singular point is 
guaranteed by the Brouwer fixed-point theorem for the 
three-dimensional Euclidean space. 

We also can prove an existence of such point from the 
electrostatic analogue. We have some bounded distri- 
bution of the positive charge density. The potential of 
the electric field is a continuous function and it decreases 
outside the charge location, tending to zero at the infin- 
ity. So, it must be maximal at some point due to the 
Bolzano Weierstrass theorem. At this point its gradient 
vanishes and the field intensity is equal to zero. This is 
a singular point. Moreover we can use the Earnshaw's 
theorem to prove that there is only two possibilities for 
the behavior of the field intensity in the vicinity of any 
singular point. Either both field intensity and the charge 
density are equal to zero inside some area, or there is 
some charge density in the singular point. We will show 
below that the later possibility leads to some dynamical 
antidrag force. 

The location of the singular point in the velocity space 
corresponds to a certain velocity. We will refer to the in- 
ertial reference frame moving with this velocity relative 
to the initial one as a "special" frame for the anisotropic 
case. Its location in the velocity space coincides with the 
singular point. Therefore, the force acting on the immo- 
bile body in this special frame due to the gravitational 
interaction with dark matter vanishes. 

In the isotropic case with nonzero u m i n the rest frame 
is surrounded in the velocity space by a set of frames in 
which the resultant force vanishes. All of them could be 
called special frames according to the definition above. 
Let us verify if in the anisotropic case there is such an 
area of special frames. Using the electrostatic analogy, 
we can reformulate this problem as follows. 

There is some static distribution of the electric charge 
with a cavity inside. Is it possible to have no electric field 
in the cavity in the case of an anisotropic charge distribu- 
tion? To answer this question we use the transformation 
of inversion r = R 2 jr. It is known that if some function 
U(r) is a solution of the Laplace equation AU = 0, then 
the function U'(r') = Rr'~ 1 U(R 2 r' /r'~ 2 ) is also its solu- 
tion [l2| • We can add an arbitrary constant to the poten- 
tial. Let us choose the value of this constant in such a way 
that we have U = inside the cavity and U = Uoo < 
at the infinity. If we choose an arbitrary point inside the 
cavity as the origin of the coordinate system and perform 
the inversion, we would turn the charge distribution in- 
side out. Note that the point charge Q — RUoo < 
appears at the point r' = 0. As a result, we get some 
spatially limited charge distribution surrounded by the 
space without electric field. This means that all the 
terms in the multipole decomposition of the field van- 
ish far from the system of charges. It is possible only 
in the case if the distribution has spherical symmetry or 
consists of spherically-symmetric parts. Thus the initial 



charge distribution before the inversion must also con- 
sist from spherically-symmetric parts. The whole charge 
distribution is not necessary spherically symmetric. An 
example of such distribution is two uniformly charged 
shells one inside another with different centers. Inside 
the inner shell there is no electric field. 

Thus an area of the special frames can exist in princi- 
ple. The distributions of dark matter velocities providing 
an area of the special frames in the velocity space can be 
thought up. But such unusual distributions can hardly 
be found in the real Universe. Indeed, we can artificially 
arrange electrical charge in a special way, but we can- 
not provide a very special distribution of velocities. One 
can hardly expect to deal with even a distribution with 
nonzero u m i n . For this reason we come to the conclusion 
that the special frame in the anisotropic case is isolated 
in the velocity space. 

If the body moves relative to the "special" frame, the 
force due to the gravitational interaction with dark mat- 
ter becomes nonzero. Let us consider the case of a small 
velocity of the body relative to the "special" frame. The 
analogy with the electrostatics is useful also for this prob- 
lem. It is generally known that the electrostatic poten- 
tial U can be expanded into a Taylor series near an arbi- 
trary point used as the origin of the Cartesian coordinates 
x,y, z 

3 3 3 

U (x, y,z) = U Q + *^2 a i x i + X! X! h ^XiX k H . (18) 

i—1 i—1 i—1 

Here Xi runs through coordinates x,y,z, Uq = £/(0,0, 0) 
is the potential at the point, and vector a is the electric 
field intensity E — —VU at the point, i.e at the origin of 
the coordinates, taken with the opposite sign. 

If the point is singular, the term with components of 
a vanishes. If we consider the potential distribution near 
the singular point we have to omit this term. The main 
coordinate-dependent term becomes the third one con- 
taining a symmetric matrix 6, which can be reduced to 
a diagonal form b = diag(6i, 62, 63) by the rotation of 
the coordinate axes. From the Poisson equation, ex- 
pressed in the Gaussian units, AU = —4irp qi we get 
61 + 62 + &3 = —2Trp q < 0. So, among these values there 
is at least one negative. 

Coming back to the considered problem we arrive to 
the conclusion that in the case when the moving frame is 
located near the special frame in the velocity space and 
its velocity V relative to the special frame is small, the 
components of the anisotropic force acting on the body 
are proportional to the components of V. If we direct 
the axes of the Cartesian coordinates x,y,z along the 
eigenvectors of the matrix b, which is a velocity-space 
analog of the matrix b appearing in (|18[) . we obtain 

F cx {hV x + b 2 V y + b 3 V z )M +b 2 + b 3 >0. (19) 

This force is similar to an anisotropic drag. In addition, 
the viscosity is negative along at least one eigenvector. 



7 



The mean friction averaged over all direction is nonposi- 
tive. 

Let us make the problem more specific. If the body 
(e.g. the Sun) rotates around the galactic center, there 
are some preferential directions, e.g. the direction to- 
wards this center, the one perpendicular to the galactic 
plane, and the one along the direction of the rotation. 
The reference frame associated with the body can be 
special first of all if the density of the particles mov- 
ing in some direction is equal to the density of particles 
moving in the opposite direction with the same velocity, 
B u (9, cj>, u) = B u (ir — 9, <f> + tt, u). This means that dark 
matter's motion is invariant relative to the inversion of 
time (T invariant). This motion can be T invariant if all 
relaxation processes are finished and if there is no rota- 
tion of dark matter relative to the galactic rotation. We 
do not know for sure anything concerning the dark mat- 
ter motion. If dark matter is rotating around the galactic 
center in different way as an ordinary matter does, the 
motions of stars are affected by small accelerations due 
to the considered effect. As a result, the ordinary mat- 
ter and dark matter subsystems of galaxy can transfer 
their angular momentum from one to another. But this 
effect is small and practically does not affect the galactic 
dynamics. 

If the velocity of the special frame relative to the frame 
of the CMB radiation isotropy is coherent with the rota- 
tion of the part of galaxy for which this frame is special, 
then the special frames for different galaxies or different 
parts of galaxy are rotating relative to each other. In 
this case there is no inertial frame, which is special for 
the whole Universe. If such a universal inertial special 
frame exists, it means that the dark matter subsystems 
of all galaxies do not rotate around their centers and the 
angular momenta of these subsystems are equal to zero. 



VI. CAPTURE OF DARK MATTER PARTICLES BY 
BLACK HOLES 



The estimation made in Sec, 3 shows that the consid- 
ered effect is significant only for supermassive black holes. 
As mentioned earlier, there is another effect acting in this 
case. Previously, we considered the dark matter without 
weak interaction to avoid any nongravitational interac- 
tion with the body. All particles penetrating the body 
escape it, as it is shown in Fig. [1] But particles penetrat- 
ing the black hole horizon cannot escape. They transfer 
their momentum to the black hole providing a force act- 
ing on it. We need a formula for the black hole capture 
cross section er for the cold dark matter to calculate this 
force. This problem was considered many times (see, e.g. 
[l3j . problem 15.11) and its solution is 



It is easy to calculate the corresponding force 



dF 



16ttG 2 M z 



dpDMU. 



(21) 



Comparing this equation with ([9]) we see that the force 
(|2"TT) is much less that (0), so we can neglect it. 

From (|20l) one can see that a minimal impact param- 
eter pmax for avoiding the capture of the dark matter 
particles is proportional to the Schwarzschild radius r g 
of the black hole p m %n — ^ r g^- For our estimation of 
the typical value of dark matter particles velocity we 
have p m i n ~ 3 x 10 3 r g . The deviation angle ip <C 1 if 

P ^ Pi = r g ) ■ If we assume that all dark matter 
particles with the impact parameter p max < p\ transfer 
their momentum to the black hole, this provide the force 



dF = 



4ttG 2 M : 



dpDMU 



(22) 



WttG 2 M 2 



-, u <C c. 



(20) 



acting on the black hole. Note that it is less that the 
force ([5]). The real force acting on the black hole due to 
the gravitational interaction with the particles with the 
impact parameters p m in < P < Pi is less that (|2"2"]l . Thus 
the lion share of the force acting on the supermassive 
black hole provide particles moving on the distances more 
that pi, where the assumption ip <C 1 holds. Also holds 
the weak gravitational field approximation and we can 
use the Newtonian mechanics. 

The major part of the force is provided by the par- 
ticles missing the body or the black hole on distances 
from about 10000 Schwarzschild radii to an R m ax many 
times more. In this region around the black hole in the 
real Universe could be an accretion disk or some orbit- 
ing compact objects. The gravitational interaction with 
them can interfere with the gravitational attraction to 
the central massive black hole, but the difference in scales 
of the effects probably cannot allow the essential devia- 
tion from the formula ([5]). The ordinary matter particles 
also can contribute to this effect. But they can inter- 
act additionally with the matter and the electromagnetic 
field surrounding the body. Ordinary matter cannot pass 
through the accretion disc around the black hole without 
nongravitational interaction, but the mirror dark matter 
can. However, the interaction of the ordinary matter ap- 
proaching the body body is much harder to study that 
the practically pure gravitational interaction of the dark 
matter particles. 



VII. CONCLUSION 

We considered the force acting on a body or a black 
hole due to the bending of dark matter particles' tra- 
jectories in the gravitational field of this massive object. 
This force depends on the velocity of the object. The 
force is also proportional to the body's mass squared, 
and the acceleration caused by it is proportional to the 
mass of the body. The estimation of its value shows that 
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it is very small even for stars. But there are objects for 
which we can expect some observational manifestation of 
the considered effect, especially on the cosmological time 
scale. These are the supermassive black holes in galaxies, 
including active galactic nuclei. But the considered effect 
can be masked by the interaction of the ordinary matter 
with other baryonic matter or the electromagnetic field 
surrounding the black hole. 

Moreover, it can be observed only in the case of es- 
sential anisotropy in the distribution of velocities of dark 
matter particles in the frame of the black hole. This effect 
can manifest in a large velocity of a supermassive black 
hole relative to other objects in the vicinity. We cannot 
measure the tangential components of the velocity and 
the effect can show itself as a difference in redshifts of 
a black hole and a part of galaxy near it. If such a dif- 
ference was detected for supermassive black holes with 
masses greater than 10 9 Mg„ n only, this would be an ar- 
gument in favor of the existence of the considered effect. 

This force allows us to select a preferable reference 
frame. In any point of the space there is a special in- 
ertial frame of reference in which there is no resultant 
force acting on the immobile body due to gravitational 
interaction with the dark matter particles flying near and 
through the body. If the distribution of the dark matter 
velocity is anisotropic in this frame, then in the frames 
moving steady and straight relative to this special frame 
there are the small forces acting on the bodies, which 
are immobile in these frames. In other words, the force 
is acting on the moving bodies in the special frame. For 
the small body velocities in the special frame, the force is 
proportional to the velocity components in the way sim- 
ilar to the case of an anisotropic drag force (fT§)) . There 
is always a directions of motion with a negative viscosity 
coefficient. The mean friction averaged over all direction 
is nonpositive. 

If the distribution of the dark matter velocity is 
isotropic in some special frame and the minimal veloc- 
ity of the particles approaching far from the body u m i n 



is equal to zero, the force acting on the body moving 
relative to this special frame with the velocity v is di- 
rected towards v. This case is similar to an isotropic 
drag with a negative viscosity. If u m in 7*= 0, then there 
are many special frames. The preferable one is the rest 
system in which the distribution of the dark matter ve- 
locities is isotropic. All the frames with the velocities 
V < u m i n relative to the rest frame are special ones. In 
other words, the is no force acting on bodies moving with 
velocities v < u min in the rest frame. If the body's ve- 
locity exceeds u m i n , the force acts towards the direction 
of motion. 

Nowadays we know some preferable inertial frames, 
e.g. the frame of CMB isotropy. To a certain degree we 
return to Aristotle's time. At any location in the Uni- 
verse there are inertial frames particular for considered 
effect. This is the special frame in the case of anisotropy 
of dark matter distribution or the rest frame for the 
isotropic case. These frames are local ones. If there is 
a universal inertial frame, which is special for the whole 
Universe, it means that the angular momenta of dark 
matter subsystems of all galaxies are equal to zero. 
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